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Introduction 


The  propagation  of  regional  seismic  waves  has  been  intensively  studied 
in  different  parts  of  the  world,  especially  in  North  America  and  Eurasia. 

Taking  the  advantage  of  a  dense  and  homogeneous  seismic  network  set  up 
in  France,  it  is  possible  to  study  in  detail  that  propagation  of  crustal  and  sub- 
crustal  waves,  by  selecting  different  local  earthquakes  well  recorded  on  this 
network. 

Attenuation  versus  distance  is  particularly  studied  and  results  obtained 
on  guided  waves  Lg  compared  with  results  on  U  S  and  USSR  structures  of 
propagation.  Q  factors  are  also  computed  versus  frequency. 

On  the  other  hand,  an  attempt  to  interprete  that  attenuation  in  terms  of 
tectonics  differences  is  proposed. 


Seismic  Data  -  Earthquakes 

25  short-period  stations  (  Z  component)  belonging  to  the  L.  D.  G.  seismic 
network  have  been  selected  for  this  study  (Figure  )  ).  The  data  are  recorded 
on  magnetic  tape  at  a  rate  of  50  samples  per  second  and  magnifications  are  gene¬ 
rally  125k  and  350k  at  1  Hz. 

The  seismometer  has  a  natural  frequency  of  1  Hz,  and  the  frequency  response  of 
the  system  ranges  from  0.5  Hz  to  20  Hz  (Figure  2). 

Earthquakes  have  been  selected  in  order  to  study  the  propagation  of  local 
waves  through  and  below  the  crust  in  France,  along  different  azimuths  of  recep¬ 
tion.  We  have  choosen  six  earthquakes  with  a  local  magnitude  around  =  4 
(  Table  I  •;  list  of  earthquakes)  located  in  or  around  France. 

Localizations  made  with  the  L.  D.  G.  network  are  given  with  an  accuracy  of 
5km  or  less  for  events  inside  France  and  up  to  10  to  20km  for  the  other  ones. 
When  possible  depths  have  been  computed.  They  are  all  shallow  but  only  three 
of  them  (West  Brittany,  Mulhouse  and  St  Pourgain  sur  Sioule)  could  be  evaluated 
with  a  reasonable  accuracy  (t  5km) 


•  / 
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Azimuths  of  propagation  from  epicenters  to  stations  are  shown  in  Figure  3. 
The  location  of  the  earthquakes  is  as  follows  : 


W.  B. 

O 

B 

U 

M 


West  Brittany 

Oviedo 

Barcelona 

Ugine 

Mulhouse 


on  the  continental  plateau 
in  the  North-West  part  of  Spain 
in  Spain 
in  the  Alps 

in  the  South  part  of  the  Rhinegraben 


Figure  4  gives  a  representative  example  of  recorded  data  and  anomalies  between 
different  travel  paths.  The  event  considered  is  the  Barcelona  earthquake  of 
12.  14.  1980  (Ml  =  3.9). 

Stations  LFF  and  LMR  which  are  located  at  the  same  distance  from  the 
epicenter  A  =  370km,  record  roughly  the  same  Pn  and  Sn  amplitudes  (LMR 
magnification  is  4  times  the  one  of  LFF).  On  the  contrary,  there  is  no  apparent 
Pg  an  Lg  waves  recorded  at  the  LMR  station,  even  though,  large  amplitude  Pg 
and  Lg  waves  are  present  at  LFF.  An  explanation  for  this  anomaly  is  that  the 
path  from  the  epicenter  to  LMR  is  oceanic  while  the  path  to  LFF  is  continental. 
The  TCF  station  located  500km  away  from  the  epicenter  receives  successively 
Pn,  Pg,  Sn  and  Lg  with  amplitudes  about  one  fourth  those  at  LFF. 


Attenuation  of  Crustal  and  Subcrustal  Phases 


a)  Global  attenuation  coefficient 

Let  us  define  A  as  the  maximum  amplitude  of  the  considered  phase  (broad 
band)  within  the  first  few  seconds  after  the  onset  time 
D  :  distance  epicenter  -  station 

^  :  global  attenuation  coefficient 

Aq  :  constant 

A  =  A0  D  *  ^ 


A 
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By  measuring  A  for  each  station,  at  Dkm  from  the  source,  it  is 
possible  to  compute  X  by  least  squares  method. 

Applied  on  the  data  described  above  the  principal  results  are  the 
following  : 

-  the  Pn  phases  is  not  fairly  well  observed  below  500  km  in  most  of 
the  cases  ;  its  global  attenuation  coefficient  %  is  of  the  order 

of  2  and  the  rather  small  dispersion  accompanying  this  result  means 
that  Pn  should  not  be  too  much  affected  by  crustal  anomalies. 

-  a  similar  attenuation  figure  is  found  out  for  the  Sn  phase,  but  as 

it  arrives  in  the  P  coda,  its  presence  is  not  always  clearly  detected. 

-  Pg  and  Lg  phases  when  present  are  strongly  energetic  in  the  range 

of  distances  of  propagation  involved  :  1  to  10°.  The  Barcelona 

earthquake  has  clearly  shown  a  lack  of  these  phases  for  oceanic 
propagation  path.  Other  data  show  that  Lg  and  Pg  do  not  propa¬ 
gate  along  paths  between  France  and  Corsica. 

Five  earthquakes  among  the  six  previously  listed  are  giving  sufficient  data 
to  draw  amplitudes  variations  versus  distance,  after  a  magnitude  normali¬ 
zation  at  Ml  =  4.  0. 

A  global  attenuation  factor  ^  around  2.  3  is  given  for  Pg  (Figure  5)  and 
a  global  attenuation  factor  Y  around  2.  5  for  Lg  (Figure  6). 

At  distances  within  100  km,  the  observed  non  linearity  of  these  two  attenuation 
curves  is  due  essentially  to  the  frequency  band  pass  of  recorded  data. 

For  both  phases,  a  large  scatter  between  earthquakes  observations  involves 
important  influence  of  azimuth,  which  can  be  interpreted  in  terms  of  crustal 
anomalies.  Nevertheless,  according  to  the  rather  small  amount  of  events 
selected  here,  all  the  suspected  anomalies  have  probably  not  been  taken  into 
account.  It  has  been  observed  that  propagation  paths  through  the  Rhone 
Valley  in  France  or  the  Ivrea  Zone  in  Italy  were  associated  with  weak  Pg  and 
Lg  wave  '  s  amplitudes  at  the  receiver  (both  zones  correspond  to  a  thin  crust 
area)  although  Pn  and  Sn  waves  are  clearly  seen  even  at  large  epicentral 
distances  (West  Brittany  recorded  in  France). 


5. 


Lg/Pg  ratio  presents  large  variations,  between  1  to  10,  with  a  mean 
value  around  3  or  4  (Figure  7).  No  correlation  with  distance  is 
clearly  seen. 

These  results  are  to  be  compared  with  those  obtained  by  Pomeroy  (1978) 
which  show  that  in  the  Western  portion  of  the  USSR  Lg  amplitudes  are 
almost  equal  to  P  amplitudes  while  in  Eastern  North  America  Lg  phases 
are  approximately  5  to  10  times  greater  than  P;  or  with  results 
obtained  by  Nuttli  (1973)  who  observed  Lg/Pg  ratio  of  the  order  of  1 
to  5  in  Western  U  S  S  R. 

X  evaluations  have  been  obtained  by  measuring  the  maximum  ampli¬ 
tudes  of  each  phase  on  broad  band  seismograms.  Every  frequency 
however  does  not  attenuate  in  the  same  way  as  it  is  clearly  shown  Figure  8 
and  Figure  9  for  the  Mulhouse  earthquake  (July  16,  1980).  After  narrow 
band  filtering  the  records  from  0.  5  to  16  Hz,  one  can  notice  a  shift  of 
amplitudes  from  high  (8-16  Hz  )  to  low  frequencies  (0.  5-0.  75  Hz  ) 
when  the  epicentral  distance  increases  from  295km  (SSF  station)  to 
518km  (RJF  station).  Amplitudes  between  Figure  9  and  Figure  8  are  in 
the  ratio  of  1  to  4. 

Consequently,  it  seems  necessary  to  analyse  this  attenuation  for  each 
frequency  band.  Such  an  analysis  shows  an  increase  of  Y  with  fre¬ 
quency  as  is  shown  for  example  in  Figure  10  for  Lg  waves  generated  by 
the  West  Brittany  earthquake.  However,  V  takes  into  account  both 
attenuation  due  to  geometrical  spreading  and  inelastic  attenuation  depending 
on  frequency. 

A  separation  between  the  two  effects  is  now  necessary. 

b)  Anelastic  attenuation  for  Lg  waves 

For  a  uniform  point  source  of  elastic  waves  in  a  spherical  earth  model, 
the  amplitudes  of  dispersed  surface  waves  measured  in  the  time  domain 
are  given  by  : 

A  =  A0  (  DsinD)*^  e’KD 
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where  Aq  is  a  factor  independant  of  distance;  D  *  ^  accounts  for 
loss  of  amplitude  due  to  dispersion;  sin_<a^  D  accounts  for  loss  of 
amplitude  due  to  geometrical  spreading,  and  e  '  describes  the  loss 

of  amplitude  due  to  anelastic  absorption;  is  1  for  waves  different 

from  the  Airy  phases  (Ewing,  1957  ). 

K  coefficient  of  anelastic  attenuation  is  expressed  by  : 

k  -  2LL 

Q  V 

with  f  frequency  of  the  wave,  Q  quality  factor,  V  group  velocity. 
In  order  to  show  the  frequency  dependence  of  this  anelastic  coefficient  K, 
we  have  represented  for  one  earthquake  West  Brittany,  amplitudes  A 
(velocities)  as  a  fonction  of  epicentral  distances  in  three  consecutive 
bandpass;  4-8Hz;  2  -  4  Hz  ;  1  -  2  Hz. 

Different  attenuation  coefficients  ^4  have  been  used  to  fit  an  attenuation 
curve  to  the  observed  values.  Going  from  o4  =  0  which  corresponds 
to  the  absence  of  attenuation  due  to  geometrical  spreading  to  'k  =  0.  5  or 
1  wich  are  tore  realistic  values,  and  finally  to  linear  attenuation 
(O  =  ,  that  is  without  anelastic  attenuation),  there  is  no  apparent  and 

clear  difference  in  the  quality  of  the  fit.  All  attenuation  curves,  however, 
show  a  strong  dependency  between  slope  and  frequency:  as  the  frequency 

increases,  the  attenuation  curve  steepons  (Figure  10).  Values  of  K  for 
Lg  waves  at  a  period  of  i  second  have  been  computed  for  each  one  of  the 
6  earthquakes,  except  the  Oviedo  event  (N-W  part  of  Spain)  for  which 
no  Lg  waves  are  recorded  in  France  (Table  II). 

The  coefficients  of  anelastic  attenuation  K  are  rather  different  from  one 
earthquake  to  another. 

The  mean  value  for  K  is  : 

K  =  0. 2  t  0.  1  d°  -  1 

which  can  be  compared  with  Nuttli  1  s  results  also  obtained  for  1  second 
Lg  waves  and  concerning  regions  at  the  same  scale  :  the  Eastern  U  S 

(  00.  7  d  '  ),  the  Northern  USSR  (  0.  1  5  d  0  “  ^  )  and  the  South  of  the 

Caspian  Sea  (  0.  35  d  0  “  1  ). 

•  A 
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The  attenuation  figure  for  France  is  comparable  to  the  one  obtained  for 
the  USSR  with  a  dispersion  large  enough  to  let  K  values  reach  the 
values  obtained  for  the  Eastern  US  (low  attenuation)  and  the  Caspian  Sea 
(high  attenuation). 

Another  order  of  comparaison  can  be  the  Western  US  (California)  where 
Richter  (1935)  obtained  0.6d°  ”  *. 

Concerning  the  rather  large  difference  between  low  K  values  for  West  Brit¬ 
tany  (  0.  16  d°  "  1  ),  Ugine  ( 0.  09  d 0  -  1  )f  Mulhouse  ( 0.  15  d°  "  1  )  and  high  K 
values  for  Barcelona  (0.29d°  "  1  )  and  S*  Pourgain  sur  Sioule  ( 0.  32  d#  -  1  ) 
an  attempt  of  intrepretation  can  be  proposed  in  terms  of  geotectonics. 

Seismic  waves  from  West  Brittany  in  the  Northwestern  part  of  France 
or  Mulhouse  and  Ugine  in  the  East  have  been  recorded  mainly  by  stations  in 
Normandy,  Morvan  and  Central  France  (North  of  the  Central  Massif)  and 
Vosges,  that  means  essentially  along  the  general  trends  of  the  Armoricain 
Massif. 

Generated  by  the  same  earthquake  in  West  Brittany,  waves  have  also  been 
recorded  in  Aquitain  region,  after  having  followed  the  trend  of  the  South 
Armoricain  accident. 

On  the  contrary,  Lg  waves  propagating  from  the  Barcelona  earthquake  to 
French  seismic  stations  have  followed  paths  towards  the  North  through  the 
Pyrenees  chain,  the  Massif  Central,  the  Sillon  Rhodanien,  and  consequently 
have  travelled  through  inhomogeneous  crustal  structures  characterized  by 
thickness  variations  and  deep  faulting  systems. 

The  propagation  for  the  S*  Pour^ain  sur  Sioule  earthquake  belongs  to  both 
of  these  two  types  of  propagation.  No  interpretation  is  proposed  here  for 
this  peculiar  event. 

Quality  Factor  Q 

Using  the  filtering  method  previously  described,  we  have  computed  Q  as  a 
fonction  of  frequency  for  each  earthquake  and  for  each  crustal  ( Figure  1 1 )  and 
subcrustal  phases  ( Figure  12  ). 
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Assuming  an  attenuation  factor  due  to  geometrical  spreading  of  0.  5,  i.  e.  : 

A  D  "  ®  e-KD 

with  :  Ilf 

K  —  - 

Q  V 

has  been  computed  as  a  fonction  of  frequency  by  least  squares  method,  the  velo¬ 
city  V  being  taken  respectively  as  : 


for 

Pn 

8.  1  km/ 

Sn 

4.  7  - 

Pg 

6.0  - 

Lg 

3.  5  - 

For  each  one  of  the  four  phases,  Q  increases  with  frequency  from  about  100  at 
0.  5  Hz  to  about  1  000  at  16  Hz,  The  dispersion  for  Pn  and  Sn  waves  is  large 
in  comparison  to  Pg  and  Lg  waves.  This  anomaly,  at  least  for  the  Q  factor 
corresponding  to  Sn  waves,  could  be  due  to  poor  Sn  amplitude  measurements  in 
a  sometimes  large  Pg  coda. 

For  Pg  and  Lg  waves,  Q  values  show  a  general  trend  of  converging  for  high 
frequencies  for  which  amplitudes  are  less  affected  by  crustal  propagation  anomalies 
than  for  low  frequencies. 

As  was  shown,  the  attenuation  coefficient  ^  due  to  geometrical  spreading  could  be 
taken  within  a  large  range  of  values  without  affecting  the  fit  to  the  data.  Each 
value  would  result  in  one  particular  Q  factor.  Consequently,  it  is  not  realistic  to 
compute  absolute  Q  values  but  only  a  general  trend  can  be  given  for  each  phase. 

The  mean  slope  of  the  Q  versus  frequency  curves  for  Pg  waves  is  around 
0.  8  -0.  9  and  represents  the  least  dispersed  figure  among  the  four  seismic  phases. 
Aki  (1975)  finds  slopes  of  0.  6  -0.  8  on  coda  waves  from  Japanese  earthquakes  and 
Nuttli  (1981)  obtains  values  around  0.  7  for  Lg  in  the  Novosibirsk  region  (USSR) 
Mitchell  (1980)  proposes  0.3-0.  5  for  Lg  waves  between  periods  of  lsec  to  40  sec. 

The  Q  versus  frequency  curves  for  Sn  waves  fall  into  the  two  family  groups  : 
one  family  (W;  Brittany  and  Oviedo)  with  'ow  attenuation  or  high  Q  values; 
one  family  (Mulhouse,  Barcelo»  id  X"  .ne )  with  higher  attenuation  or  lower 
Q  values.  In  the  first  case  clear  Sn  data  have  been  recorded;  on  the  contrary 
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weak  and  doubtfull  Sr  arrivals  possibly  mixed  up  with  Pg  wave  coda  are  seen  in 
the  second  case. 

As  was  printed  out  formerly,  the  attenuation  for  the  West  Brittany  earthquake 
appears  to  be  the  lowest  one  and  gives  high  Q  values  for  each  of  the  four  phases. 


Conclusion 


Using  the  seismic  data  from  one  homogeneous  short  period  network  installed 
in  France,  we  have  studied  the  propagation  of  regional  crustal  phases. 

Global  attenuation  factors  on  broad  band  data,  taking  into  account  both  geometrical 
spreading  and  anelastic  attenuation  are  of  the  order  of  2,  for  Pn  and  Sn  waves. 
They  are  sligthly  higher  for  Pg  (  2.  3  )  and  Lg  (2.  5)  waves  but  the  amplitude  of 
these  waves  is  largely  dependent  on  the  azimuth  0i  propagation.  It  has  been  ob¬ 
served  that  Pg  and  Lg  are  strongly  attenuated  when  propagating  through  the 
Sillon  Rhodanien  in  France  or  the  Ivrea  zone  in  Italy,  although  Pn  and  Sn  phases 
are  clearly  received  for  the  same  paths. 

Lg  /  Pg  ratios  present,  as  a  consequence,  a  large  scatter  between  1  and 
JO,  with  a  mean  trend  around  3  to  4. 

-/  using  filtering  techniques  we  were  able  to  pursue  the  study  further  and 
evaluate  anelastic  attenuation  and  quality  factor  for  different  frequencies. 

We  found  that  the  Lg  attenuation  factor  is  of  the  order  of  0.  2  d'  '  1  at  1  Hz, 
a  result  to  be  compared  with  the  values  of  0.  07  d  0  "  ^  ( Eastern  US),  0.  1  5  d  °  "  * 
and  0.  35  d°  *  ■  (  respectively  Northern  USSR  and  South  of  Caspian  Sea)  obtained 

by  Nuttli.  Apparently,  the  attenuation  values  for  France  andNorthern  USSR 
are  similar. 

Q  factors,  computed  for  each  of  the  local  phases,  show  a  clear  dependence  on 
frequency,  even  if  their  absolute  value  cannot  be  confidently  assessed. 

One  must  nevertheless  be  carefull  about  interpretating  this  frequency  depen¬ 
dence  of  Q.  Laboratory  studies  on  rock  samples  have  not  pointed  out  any  vari¬ 
ation  of  Q  factor,  except  very  recently,  see  Bonner  and  al098J). 

The  apparent  variations  of  Q  obtained  in  this  paper,  might  be  a  result  of 
scattering  effects  or  of  different  propagation  paths  associated  with  different 
frequencies. 
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FIGURE  CAPTIONS 


Figure 


Figure 


Figure 


_1_  25  short  period  seismic  stations  (vertical  component)  which  data 

are  telemetered  to  a  central  station  close  to  Paris  where  they  are  recorded 
on  magtape  (digital).  Data  are  sampled  at  50  samples/sec,  with  magni¬ 
fications  of  125  000  to  350  000  at  1  Hz. 

The  seismometer  has  a  proper  frequency  of  1  Hz.  The  frequency  band 
recorded  on  magnetic  tape  is  0.  5  -  20  Hz. 

2^  Seismometer  frequency  response  in  terms  of  magnification 
©:  seismometer  alone 

©:  seismometer  with  amplifier  and  paper  recorder 
_3  Propagation  paths  from  six  earthquakes  to  French  seismic  stations 
W  B  :  West  Brittany  on  the  continental  plateau 
O  :  Oviedo  in  North  West  part  of  Spain 
B  :  Barcelona  in  North  East  part  of  Spain 
U  :  Ugine  in  the  northern  Alps 
M  :  Mulhouse  in  the  southern  Rheingraben 


Figure 


Broad  band  records  in  3  seismic  stations 


LFF  and  LMR 


(A  =  370km);  TCF  (  A  =  500km)  for  the  Barcelona  earthquake, 
magnifications  (at  1  Hz )  are  32  000  (LFF)  and  125  000  (LMR,  TCF). 

Tick  marks  every  second.  Pg  and  Lg  phases  are  energetic  when  recorded 
at  LFF  (  continental  path)  ;  but  no  Pg  and  weak  Lg  are  recorded  at  LMR 
at  the  same  distance  (oceanic  path). 

Pn  and  Sn  are  clear  at  both  stations  LFF  and  LMR  ;  very  weak  at  TCF. 

[  Peak  to  Peak  amplitudes  of  Pg  waves  versus  distance.  2.  3  is  the 
global  attenuation  coefficient.  Data  are  taken  from  broad  band  seismograms, 
i  Same  figure  caption  for  Lg  phases. 

Lg/Pg  ratio  versus  distance. 

Figure  8  and  Figure  9  Digital  filtering  process  is  applied  by  using  FIR  technic 
(  Finite  Impulse  Response)  for  different  frequency  bands.  At  the  bottom 
of  each  figure  is  the  broad  band  signal. 


Figure 


Figure 

Figure 
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Figure  10  For  three  band  passes,  attenuation  curves  with  different  geometrical 
spreading  coefficients  o,  ,  are  adjusted  on  data. 

Figure  11  and  Figure  12  For  Pn,  Sn  and  Pg,  Sg  respectively  quality  factors 
are  computed  versus  frequency  for  six  earthquakes  (when  these  phases 
exist ). 

A  coefficient  of  geometrical  spreading  equal  to  =  0.  5  has  been  taken 

for  that  purpose. 


A 


<n 


Table 


INELASTIC  ATTENUATION  FOR  Lg  (I Hz) 


A  =  Ao#(D  sinDXl/2*EXP(-KD) 


EARTHQUAKES 

K  Cd 

W-BRITANNY 

0.16 

OVIEDO 

no  Lg 

BARCELONA 

0.29 

UGINE 

0.09 

MULHOUSE 

0.15 

St  POURCAIN 

0.32 

Table  II 


980 


o 


QUAKE  OF  MULHOUSE  -  JULY  16,1980 


Figure 


QUAKE  OF  MULHOUSE  -  JUL 


Figure 


Lg  ATTENUATION  VERSUS  DISTANCE  FOR  WEST  BRITTANY  EARTHQUAKE 


•Wa 


I 
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Abstract 


"^Global  attenuation  versus  distance  has  been  determined  for  each  phase  Pn, 
Sn,  Pg  and  Lg  and  the  influence  of  frequency  analysed.  For  Pn  and  Sn 
waves  they  are  of  the  order  of  2,  and  slightly  higher  for  Pg  (2.  3)  and  Lg 
(2.  5)  waves.  But  the  amplitude  of  these  waves  is  largely  dependent  on  the 
azimuth  of  propagation.  . 

Lg/Pg  ratios  present  a  large  scatter  between  1  and  10  with  a  mean  trend 
around  3  to  4.  A.IO.l ' 

We  found  that  the  Lg  anelastic  factor  is  of  tne  order  of  0.-2'^“-“'*  at  1  Hz, 
a  result  to  be  compared  with  the  values  of  0.-0^-** — &  (Eastern  US), 
0^T5d*“  *  and  (respectively  Northern  U  S SR  and  South  Caspian 

Sea),  obtainoif  by  NuttlM ) . 

£^The  < quality  factor,  computed  for  each  of  these  local  phases,  shows  a  clear 
dependence  on  frequency,  increasing  from  some  100  at  1  Hz  to  some  1  000 
at  10  Hz  probably  due  to  scattering  effects. 
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